% 1. stance knee
load('/home/shu/workspace/Research/2DKnee_Torso_Outputs/data/Afit1.mat')
ns_time = mean(time_rec)*ave_timeSK';
% [x_meanSK,upperBSK,lowerBSK]=meanValue(ave_dataSK_s);
meanValuePlotG1(ave_timeSK,upperBSK,lowerBSK,x_meanSK,'SKnee.eps',...
    'Angle(rad)',a0(3,:),2,ns_time,6); 
%%
% 2. non-stance knee
[x_meanNSK,upperBNSK,lowerBNSK]=meanValue(ave_dataNSK_s);
meanValuePlotG1(ave_timeNSK,upperBNSK,lowerBNSK,x_meanNSK,'NSKnee.eps',...
    'Angle(rad)',a0(4,:),2,ns_time,7)
% meanValuePlot(ave_timeNSK,upperBNSK,lowerBNSK,x_meanNSK);


% 3. non stance slope
[x_meanNSL,upperBNSL,lowerBNSL]=meanValue(ave_dataNSL_s);
% meanValuePlot(ave_timeNSL,upperBNSL,lowerBNSL,x_meanNSL);
meanValuePlotG1(ave_timeNSL,upperBNSL,lowerBNSL,x_meanNSL,'NSSlope.eps',...
    'Slope',a0(2,:),2,ns_time,3)


% 7. hip position
[x_meanHP,upperBHP,lowerBHP]=meanValue(ave_dataHP_s);
meanValuePlotG1(ave_timeHP,upperBHP,lowerBHP,x_meanHP,'HipPos.eps',...
    'Position(m)',a0(1,:),1,ns_time,1);